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Abstract. As all the globally fitted nuclear PDFs (nPDFs) have been so far impact parameter 
independent, it has not been possible to calculate the hard process cross sections in different 
centrality classes consistently with the global analyses. In [T] we have offered a solution to 
this problem by determining two spatially dependent nPDF sets, EPS09s and EKS98s, using 
the ^-systematics of the earlier global fits EPS09 and EKS98 and an assumption that the 
spatial dependence can be written as a power series of the nuclear thickness function. For a 
data comparison, we have calculated the nuclear modification factor of inclusive neutral pion 
production in d+Au collisions at RHIC in four centrality bins at midrapidity and compared 
these to a PHENIX measurement. In addition, we have also performed a similar calculation for 
inclusive photon production in d+Au collisions at RHIC. 



1. Introduction 

In the collinear factorization framework [21 [3] the cross section for inclusive k production in a 
hard process of a heavy ion collision between nuclei A and B can be computed as 

da AB^k+x = J- tf(x,Q 2 )® ff{x,Q 2 )®da^ k+x ' + 0(1/Q 2 ), (1) 

i,j,X> 

where f± {ff) is the process-independent nuclear parton distribution function (nPDF) for a 
parton flavor i (j) in the nucleus A (B) and da i 3~ yk ~*~ x represents the partonic pieces which 
can be computed using perturbative QCD (pQCD). Similarly as the free proton PDFs, also 
the nPDFs can be determined via a global analysis considering as many different processes as 
possible. Usually the nPDFs are determined in terms of the free proton PDFs, , (x, Q 2 ) and 
the nuclear modification of the PDFs, Rf(x,Q 2 ), as 

f t A (x,Q 2 ) = Rf(x,Q 2 )f l N (x,Q 2 ). (2) 

So far all the globally analyzed nPDFs have been spatially independent. However, it is reasonable 
to assume that as the nucleus itself is not spatially uniform, also the nuclear modifications of 
the PDFs should somehow depend on the position of the nucleoli inside the nucleus. In Ref. [1] 
we have addressed this issue by considering the A-systematics of two globally fitted nPDF sets, 
EPS09 (3j and EKS98 [5j using the framework discussed in the next section. 



2. Framework 

First we introduce a nuclear modification of the PDFs, rf(x, Q 2 ,s), which now depends also on 
the transverse position s of the nucleon inside the nucleus. For this we require that the spatial 
average of the quantity gives back the spatially independent modification, 

R?(x,Q 2 ) = j J d 2 sT A (s)rf(x,Q 2 ,s), (3) 

for which we take the values from earlier global fits (EKS98 or EPS09). This, however, does 
not restrict the form of the spatial dependence in any way, so we have to make an assumption 
for that. Motivated by the shadowing region at small x, we assume that the rf(x, Q 2 ,s) can be 
written as a power series of the nuclear thickness function T^(s): 

n 

rf(x, Q 2 , s) = 1 + £ ^(x, Q 2 ) [T A (s)f , (4) 
j=i 

where now the fit parameters dj(x,Q 2 ) for each parton flavor depend only on x and the scale 
Q 2 but not on A. This is important for correct mapping of the spatial dependence with the 
^-dependence of the globally fitted Rf. In practice we obtain the values for our fit parameters 
dj(x, Q 2 ) by minimizing the \ 2 defined as 

Rf(x,Q 2 ) - \ Jd 2 sT A ( S )rf(x,Q 2 ,s) ] 2 

Wfix^Q 2 ) \ K ' 

in a (x, Q 2 ) grid for all the parton flavors. The weight factor W^(x, Q 2 ) is set to 1 (1 — R.f(x, Q 2 )) 
for the EPS09 (EKS98) analysis. For the EPS09 analysis we perform this fitting also for the 30 
error sets both in leading (LO) and next-to-leading order (NLO). As can be seen from figure [IJ 
the A-systematics of Rf is very well reproduced with the power series ansatz when we consider 
A > 16. In our analysis we found out that the first four non-trivial terms of the power series in 
equation |4] are enough for an accurate fitting in the whole x and Q 2 region, for all the different 
sets considered. 

After the values for the fit parameters c* (x, Q 2 ) are obtained through this fitting procedure, 
we can calculate the rf(x, Q 2 ,s) and determine the new spatially dependent nuclear modification 
sets, EPS09s and EKS98s, in which "s" stands for "spatial". These are now published in our 
websit^] To illustrate the spatial dependence of the new nPDFs, the gluon modification from 
EPS09sNLO for lead nucleus is plotted in figure [2] as a function of x and the transverse distance 
s = |s|. The general feature is that the nuclear effects are larger in the dense center (small s) of 
the nucleus and smaller at the sparse edge (large s). 



xHx,q 2 ) = Y. 

A 



3. Applications 

The nuclear effects of an observable can be quantified using the nuclear modification factor 
Raa- The improvement here is that now with the new nPDF sets one can compute the Raa 
also for different centrality classes in a manner which is consistent with the globally fitted nuclear 
modifications of the PDFs. Detailed instructions of the implementation of these new nPDFs 
can be found in Ref. pQ. 

In figure 3] we have plotted the nuclear modification factor R^Am f° r inclusive 7r° production 
in d+Au collisions at ^/snn = 200 GeV and y = as a function of px in four centrality classes. 
The calculations are done in NLO (with the INCNLO-packagfQ using the CTEQ6M PDFs [6] 
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Figure 1. The ^4-dependence of Rf(x,Q 2 ) at 
fixed x and Q 2 values from the central sets of 
EPS09NLO (crosses) and LO (pluses) and from 
EPS09sNLO (green) and LO (blue). From pp. 



Figure 2. The gluon modification in a lead 
nucleus r^ h (x,Q 2 ,s) from EPS09sNLO as a 
function of x and s at the EPS09 initial scale. 
From pp. 



with the EPS09s modifications pp and three different fragmentation functions (FFs), KKP [7], 
AKK [8j, and fDSS (9]. The uncertainty bands are calculated from the error sets of EPS09s with 
the fDSS FFs. The centrality classes are defined in terms of impact parameter intervals, which 
are calculated using the optical Glauber model. The PHENIX datapoints [10J have been scaled 
by overall factors which all are consistent with the overall normalization uncertainties given by 
the experiment. The corresponding results for the forthcoming p+Pb collisions at the LHC can 
be found in pQ. 
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Figure 3. The nuclear modification factor for inclusive ir° production in d+Au collisions at 
RHIC in four centrality classes at midrapidity, calculated with the EPS09s nPDFs and different 
FFs. The data are from PHENIX [TO]. From pp. 



We have also performed similar calculations for inclusive photon production, -R^Au' m 
d+Au collisions at RHIC, which are shown in figure [4} The results for the minimum bias 
nuclear modification factor for this process, calculated with different nPDFs, can be found in 
Ref. The inclusive photons consists of two components: direct and fragmentation. For the 
fragmentation component we have used the BFG (set II) FFs |12j . otherwise the setup is the 
same as for the 7r°'s above. To quantify the isospin effect, we have also plotted the FQ Au using 
the free nucleon PDFs in each panel. Similarly as for pions, also for the photons we can see that 
the nuclear modifications are larger for central collisions than for peripheral collisions. 
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Figure 4. The nuclear modification factor for inclusive 7 production in d+Au collisions at 
RHIC in four centrality classes at midrapidity, computed with nPDFs (EPS09s) and proton 
PDFs (CTEQ6M). From p]. 
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